Energy of Magnetic Vortices in Rotating Superconductor 
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Abstract. We carry out a systematic analytic investigation of stationary and cylindrically sym- 
metric vortex configurations for simple models representing an incompressible non-relativistic su- 
perconductor in a background, which is rigidly rotating with the angular velocity Jloo- It is shown 
that although the magnetic and kinetic contributions to the energy per unit length of such a vortex 
are separately modified by the background angular velocity, its efi'ect on the total energy per unit 
length cancels out. For a type II superconductor threaded by a parallel array of such vortices, this 
result implies that the relevant macroscopic magnetic field strength H will not be equal to the large 
scale average (B) of the local magnetic induction B (as has previously been suggested) but instead 
that H will simply be equal to the external London field value Boo = — (2m/q)iloo (where m and q 
axe the mass and charge of the condensate particles) that characterizes the value of B outside the 
vortices. 



I. INTRODUCTION 

Because of the property of irrotationality of superfluid- 
ity, it is well known that the angular momentum of a zero 
temperature superfluid is "carried" only by the vortices 
it contains. Similarly, when a magnetic field is applied 
to a type II superconductor, there exist quantized flux 
tubes. When one considers a rotating superconductor, 
the two effects, velocity and electromagnetism, are com- 
bined. It is the purpose of this work to carefully analyze 
the energy of such a system. 

The physical motivation behind this work is the need to 
clear up some confusion that has arisen in the context of 
neutron star mattend about the relation between the large 
scale magnetic field strength H and the average value (B) 
of the local nragnetic induction _B in a rotating type II 
superconductor threaded by a parallel array of vortices. 
In order to clarify the issue we shall proceed on the basis 
of the same kind of simplification that was postulated 
as the basis of the earlier discussion, lil working in terms 
of a broad category of non-relativistic incompressible su- 
perfluid models that includes, but is not restricted to, 
the special case characterized by the standard Ginzburg- 
Landau ansatz. The main conclusion of our work, as will 
be shown in Section [V|, is that the conventionally defined 
macroscopic field strength H will simply be given by the 
value (g) of the external (London) limit B^o of the local 
induction B, and not, by the average value (B) as was 
previously suggested. El 

The essential feature of the models to be dealt with 
is just the usual postulate that the relevant charged su- 
perfluid constituent is represented by a locally variable 
number density rig of bosonic particles that are charac- 
terized by an effective mass m, charge q and a momentum 
covector having space components 



where h is the Dirac-Planck constant, Ai is the magnetic 
vector potential and (p is a scalar with period 27r repre- 
senting the phase variable of the boson condensate. In 
ordinary laboratory applications the particles would rep- 
resent Cooper type electron pairs, characterized in terms 
of the charge and mass of the electron by the exact rela- 
tion q — — 2e, and to a good approximation by tti ~ 2me, 
whereas in the context of neutron star matter they would 
represent proton pairs, characterized hy q = 2e and an ef- 
fective mass given roughly in terms of that of the proton 
by m « 2mp. 

The scenarios we shall consider will be of the usual 
kind, in which each individual vortex is treated as a sta- 
tionary cylindrically symmetric configuration consisting 
of a rigidly rotating background nrediunr with uniform 
angular velocity floo , say, together with a charged super- 
fiuid constituent in a state of differential rotation with 
a velocity v, which tends at large distance towards the 
rigid rotation value given by rioo?", where r is the cylin- 
drical radial distance from the axis . It will be supposed 
that the superfiuid particle number density ris vanishes 
on the axis and is a monotonically increasing function 
of the cylindrical radius variable r, tending rapidly to 
a constant value rioo at large distances from the axis: 
JT-s = JT-oo for r ^ ^, say, where ^ is a parameter inter- 
pretable as the core radius. It will further be supposed 
that the local charge density is canceled by the back- 
ground so that there is no electric field, but that there 
is a magiretic induction field with magnitude B and di- 
rection parallel to the axis, whose source is the axially 
oriented electromagnetic current whose magnitude j will 
be given by 



j = qns{v - floor). 



(2) 



+ qAi = h\/i(f, 



(1) 



The relevant Maxwellian source equation for the mag- 
netic field will have the familiar form 



dB 

dr 



-— = -47rj . 



(3) 



The other relevant MaxweUian equation is the one gov- 
erning the axial component A (which in an appropriate 
gauge will be the only one) of the electromagnetic po- 
tential covector, which will be related to the magnetic 
induction by 



d{rA) 
dr 



rB . 



(4) 



The essential property distinguishing the "superconduct- 
ing case" from its "normal" analogue is the London 
flux quantization condition, which in the present context 
(where all physically relevant quantities depend only on 
the cylindriCjal radius r) will be expressible in the well 
known formtl 



mv + qA 



Nh 



(5) 



where N is the relevant phase winding number, which 
must be an integer. 

It is to be noted of course that by themselves the fore- 
going equations are not quite sufficient to fully charac- 
terize the model: in order to obtain a complete system 
it is also necessary to have some well defined prescrip- 
tion for the radial dependence of the number density n^ , 
which will be referred to below as the structure function. 
The available literature does not seem to provide any 
fully adequate general purpose ansatz for such a struc- 
ture function, though various, more or less satisfactory, 
phenomenological prescriptions have been put forward in 
particular contexts. One of the simplest proposals is to 
postulate that n^ falls discontinuously from its asymp- 
totic constant value noo to zero. Such a simple ansatz is 
in fact perfectly adequate for many purposes, since, as 
will be seen below, much of the relevant physics turns 
out to be insensitive to the detailed structure of the core. 
However, no such specific prescription for the structure 
function will be needed to obtain the general result of 
section ^ 

The plan of this paper is the following. In section 
P, we transform our system of equations to a simpler 
form by considering the deviations of all quantities with 
respect to their asymptotic values corresponding to rigid 
rotation. Section HI is devoted to the demonstration of 



the cancellation between the rotation-induced terms of 
the kinetic energy and the magnetic energy. In section 
[rv| , we show that there is a simple relation between the 
total energy per unit length of the vortex and the total 
flux independently of the details of the structure of the 
vortex. Section VI and VII are concerned respectively 
with the external and the internal solution representing 
the vortex. In section |VIII we evaluate explicitly the 
energy contributions as functions of the core parameters 
and finally, in section M we apply our results for one 
vortex to the case of an array of aligned vortices and 
obtain our main result concerning the macroscopic field 
strength H. Section IX summarizes this work. 



II. HOMOGENIZATION OF THE SYSTEM 

For given values of the relevant physical constants m 
and q and the rotation rate fioo, and subject to the pro- 
vision that the structure function for ris has been pre- 
scribed in advance, the foregoing equations will consti- 
tute a linear differential system relating the variable func- 
tions V, B, A to the integer valued parameter N. Before 
proceeding, it will be useful to take advantage of the pos- 
sibility of transforming the preceding system of equations 
to a form that is not just linear but also homogeneous by 
replacing the variables v, B, A by corresponding vari- 
ables V, B, A that are defined by 



V 



n.. 



B = B-Ba 



A = A- IrB^ 



(6) 



(7) 



(8) 



Here B^ is the uniform background magnetic field value 
that would be generated by a rigidly rotating supercon- 
ductor and is given by the London formula. 






(9) 



obtained by combining (p[) and (H) in the specialized case 
of rigid corotation, which is V = 0. 

In terms of these new variables the equation (Q) will 
be transformed to the form 



d{rA) 
dr 



rB, 



(10) 



while the other differential equation (g) will be trans- 
formed to the form 



dB 

dr 



Anj , 



in which we shall have 

j = qrisV . 



(11) 



(12) 



Finally the flux quantization condition (pf) will be con- 
verted to the form 



mV + qA= , 

r 

which can be used to transform (HO) to the form 
771 d{Vr) 



qr dr 



B. 



(13) 



(14) 



The advantage of this reformulation is that unlike v, B 
and A, the new variables V, B and A are subject just to 
homogeneous boundary conditions: they must all tend to 
zero as r — > oo, while at the inner boundary, as r — > 0, 



there is just the regularity requirement that B should 
be bounded, so that we have B —^ Bq for some finite 
limit value Bq, which by (10) entails automatically that 
A should tend to zero. Since the number density n^ is 
postulated to vanish at the origin there is no correspond- 
ing restriction on V. We have thus obtained a homoge- 
neous linear system of equations relating the integer N 
to the set of three functions consisting of the excess (with 
respect to the background) magnetic induction variable 
B, and the corresponding excess potential variable A to- 
gether with the relative velocity variable V, or equiva- 
lently the current magnitude j as given by (|l2). This 
means that they will be expressible in the form 



B^NB, 



V = NV, 



A^NA, 



J = Nj , 



(15) 



(16) 



in terms of corresponding rescaled functions B, A, V and 
j that will be fully determined (independently not just 
of the rotation parameter Boo = —2mfloo/q but also of 
the winding number N) just by the physical constants 
771 and q and the specification of the structure function 
giving the radial dependence of the number density Us . 



while fkin is the corresponding deviation of the kinetic 
energy from that of the state of rigid corotation at the 
angular velocity rJoo, namely 



_ m 



{v' - ^Ir^) 



(21) 



Note that in addition to f/mag and [/kin the total energy 
per unit length C/tubc associated with the vortex will con- 
tain an extra potential energy term allowing for effect of 
the breakdown of superfluid condensation in the core, 
but this will not be relevant for the work of the present 
section. In the limiting case of an ordinary superfluid, as 
characterized by vanishing charge 9 = 0, the kinetic con- 
tribution would be the dominant one, but in the context 
of sup|erconductivity, i.e., when q is non-zero, it is com- 
monlylil overlooked, perhaps because of the small value of 
the electron mass that is relevant in laboratory applica- 
tions. The purpose of the present section is to show not 
only that the kinetic contribution will not in general be 
negligible compared with the magnetic contribution, but 
also that its inclusion brings about considerable simplifi- 
cation. 

To start with, using the the decomposition (0) of the 
magnetic field, it will be possible to express the magnetic 
energy density contribution in the form 



III. ROTATION ENERGY CANCELLATION 
LEMMA 

One of the main purposes of the present work is to 
demonstrate, in the present section, a useful lemma con- 
cerning mutual cancellation - independently of the radial 
dependence of the relevant particle density Ug - between 
the background rotation dependent term in the magnetic 
energy per unit length 



t^mag — / ^raagdb , (1') 



and the corresponding term in the kinetic energy per unit 
length 



C/ki, 



C-kin "O . 



with 



dS — 2nr dr . 



(18) 



(19) 



In the above expressions, fmag is the extra magnetic en- 
ergy density arising from a non-zero value of the phase 
winding number N, i.e., the local deviation from the mag- 
netic energy density due just to the uniform field i?oo £^s- 
sociated with the state of rigid corotation at the angular 
velocity 17 oo, namely 



tmag - g^ 






(20) 



imag - g^ 



BooB 

An 



(22) 



while similarly, using the decomposition (o) of the ve- 
locity, it will be possible to express the corresponding 
kinetic energy density in the analogous form 



£kin = yr7,(v2 + 2r!ooVr), 



(23) 



which can usefully be rewritten in terms of the current 
magnitude j, using (O) and M) as 



C-kin = y"-sV 



2 - 



(24) 



Using the Maxwell source equation ( |ll| ) this can be con- 
verted to the form 



^kii 



-nsV' 



Boo d ^r^B 
r dr \ 8tt 



BooB 

An 



(25) 



in which the second term can be seen to be a pure di- 
vergence, while the last term can be seen to be equal in 
magnitude but opposite in sign to the last term in (p2), 
so that there will be a cancellation between them when 
the magnetic and kinetic contributions are combined. 

At an integrated level, in view of (|2^), it will be possi- 
ble to express the magnetic energy in terms of quantities 
[/mag and <& that are specified independently of Boo, in 
the form 



f/, 



maj:^ — ^mag 



Bo 

An 



-<!>, 



(26) 



where the part that would still be present if the back- 
ground were non rotating is given by 



C/n 



fi2 



(27) 



and where the coefficient $ is a flux integral of the simple 
form 



$= / BdS^N(t), 
where (j) is the usual flux quantum, given by 

Q 



(28) 



(29) 



In a similar manner, it will be possible to express the 
kinetic contribution in terms of quantities C/kin and <i>kin 
that are also specified independently of fioo, or equiva- 
lently of Bac , in the form 



C^kin = C^kin + -; $kin ■ 



(30) 



where the part that would still be present if the back- 
ground were non rotating is given by 



(31) 



C/kin = ^ I jVdS, 



and where the coefficient $kin is given by 

$ki„ = t: j r^'^dr = TT j d{r^B) - f 2TrrB dr , (32) 

which corresponds to the last two terms on the right hand 
side of (|25|). The first term in this expression clearly 
vanishes when the integration is taken over the whole 
range from the center, where r = 0, to the large radius 
limit where r^B — > 0, as can be seen from the explicit 
solution (pq). We are thus left with the second term, 
the kinetic analogue of the magnetic flux contribution, 
to which it is evidently equal in magnitude but opposite 
in sign, i.e., we obtain 



$kin = -$. 



(33) 



It can thus be seen that there is a remarkable cancellation 
whereby the dependence on Boo, or equivalently on fJoo, 
in the separate magnetic and kinetic energy contributions 
will cancel out when they are combined, so that we are 
left simply with a result of the form 



^mag 1 (Vki] 



t/„ 



C/ki, 



(34) 



Since the terms in this expression are both quadratically 
dependent on fields, namely V and B^ that by (03) and 



(Eq) will just be proportional to the winding number TV, 
we obtain the following conclusion. 

Rotation energy cancellation lemma: Whereas 
the separate values of the the magnetic and kinetic con- 
tributions (as defined using the formulae (pq) and (pfl) 
above) to the energy per unit length of the vortex will 
be affected by the rate of rotation of the background fioo 
(or equivalently the corresponding London field -Boo — 
—2jndao/q), the combination of these two contributions 
will not depend directly on l^oo and can be simply ex- 
pressed in the form (as a result of equations (fsl) and 



Un 



C/ki, 



UN-" 



(35) 



where we recall that N is the winding number and U 
depends only on the physical constants m and q and on 
the form of the radial distribution of the number density 
Us- A simple form for U will be given in the next section. 



IV. AXIS-FIELD ENERGY FORMULA 

The preceding result, namely the cancellation of the 
contributions due to the background rotation, was ob- 
tained simply with the background London equation (|9|) 
without using the full, i.e., local London quantization 
condition (|^) . The ultimate cancellation of the B^o - de- 
pendent contribution is attributable at a local level to 
the fact that the -Boo ~ dependent contribution to the 
combined energy density is a pure divergence: 

87r(£,„ag + £kin) =B' + 4^-jV + ^^{r'B) . (36) 
\ ° / q r dr 

We can obtain a stronger result if we now invoke the 



more specialized relation (n4 
the quantization condition 



which is a consequence of 
that specifically charac- 



terizes superconductivity. This condition can be seen to 
imply that the whole of the right hand side of ( |36| ) will 
be expressible as a divergence, since we shall have 



B^+A7:-jV = --^(rVB). 
7 or dr 



q 



(37) 



It can thereby be seen, using dq) again, that the combined 
energy density will be expressible as 



£n 



f kii 



1 d 

Snr dr 



((.. 



rA- 



'i> 



(38) 



In the outer limit, as r — > cx), the rapid fall off of B 
will ensure that the quantity inside the divergence will 
tend to zero. In the inner limit, as r ^ 0, the first term 
in the divergence obviously gives no contribution, and 
the consideration that A should be bounded ensures that 
the second term also gives no contribution, so we shall 
be left with the contribution just from the final term, 
which is proportional to the winding number N . The 



final outcome of the integration of ( p8| ) can be stated as 
follows. 

Axis-field energy lemma: Subject to the London 
quantization (as given by (|5|) above) the combination of 
the magnetic and kinetic contributions (as defined using 
the formulae (26) and (^) above) to the energy per unit 
length for a vortex with given winding number N and 
corresponding total fiux <i> as specified by (28) will be 
provided just by the the axis-field value Bq according to 
the proportionality law 



Ur, 



C/ki, 



8n 



(39) 



where Bq is the value on the axis of the relative magnetic 
field value B as given by (M), i.e. it is the difference 



Bo^Bo- B, 



(40) 



between the central value, Bq, of the magnetic induction 
B and its asymptotic London value i?oo- 

A corollary of this second lemma is that the combina- 
tion of the kinetic and magnetic energy per unit length 
will remain the same whatever the internal structure, as 
long as the total flux and the axis magnetic field are the 
same. The simplest such configuration is given by a field 
B retaining the same uniform central value Bq out to a 
cut-off where it drops discontinuously to its asymptotic 
value Boo- This cut-off radius, R, say, is adjusted so as 
to give the same total fiux as in the actual model, i.e. so 
as to satisfy the specification 



Bq Bq 



(41) 



Since the quantity Bq , i.e. the value on the axis of the 
rescaled field defined by (|l5|), depends only on the phys- 
ical constants m and q and on the form of the structure 
function specifying the radial dependence of the number 
density Ug, it follows that the same applies to the effective 
radius R, which will thus be independent of N , as well 
as of the background rotation rate fico = —qBoo/2m. 
The conclusion that the effective magnetic radius de- 
pends only on the structure function specifying n^ is in- 
terpretable as a restatement of our first lemma, since it 
can be seen that the coefficient IJ in ( p5| ) will be given 
just in terms of this effective radius R by the formula 



U = 



2q'^R' 



(42) 



AVERAGE OVER AN ARRAY OF ALIGNED 
VORTICES 



Let us now consider the typical situation in a type II 
superconductor, in which we have not just a single vortex 



but a parallel array of such vortices with sufficiently low 
mean number density per unit surface area, v, say, for 
the separation distance between neighboring vortices to 
be large compared with the penetration length A. Since, 
according to (R) and (Eq), each vortex carries an extra 
magnetic flux $ in addition to the contribution from the 
uniform London field -Boo, the large scale average mag- 
netic field will be given by 



{B) - So 



1^$. 



(43) 



As compared with the average energy density of a config- 
uration in rigid corotation with the given angular velocity 
r^oo , but with no magnetic field, the extra energy density 
averaged over a large number of vortices will be given by 



{£) — ^Lon + (^tu' 



be/ 



(44) 



where ^Lon is the uniform contribution from the London 
magnetic field, i.e. 



^Lon — 






(45) 



and where (i?tubo) is the large scale average of the contri- 
bution given locally for the separate vortices by 



i^tuh 






'-'kin ~r *-'Ci 



(46) 



where ^mag and £^kin are the magnetic and kinetic energy 
contributions discussed in the preceding sections and i?con 
is the condensation energy contribution depending just 
on the radial distribution of the condensate number den- 
sity Us (in a manner that is unimportant for our present 
purpose) , which means that the corresponding additional 
contribution 



t/c, 



^coii '^■O 



(47) 



to the vortex energy per unit length can be treated just 
as a constant as far as the present section is concerned. 
It follows that we shall have 



(^tube) = VU , 



(48) 



where U is the total energy per unit length of an individ- 
ual vortex as given by the combination 



[/=[/, 



mag 



C/ki, 



C/o 



(49) 



in which the first two contributions will separately de- 
pend on the background rotation velocity Sloo (or equiv- 
alently on the London field Boo) but in which, by the 
cancellation lemma expressed by (pi]), the total (like fi- 
nal term) will not. 

Since each vortex is associated with a momentum cir- 
culation of magnitude 2nhN there will be a correspond- 
ing generalized vorticity, in the sense of momentum cir- 
culation per unit area, with large scale average given by 



{w) = 2TrhNiy. (50) 

In terms of this quantity the large scale average ( |44| ) of 
the extra energy due to the deviation from a configura- 
tion of unmagnetized rigid corotation will be given by 



o2 



(w) U 



8tt 2TTn N 



(51) 



For a large scale variational description it is convenient to 
use (w) and (B) as the independent variables. In terms 
of these, the London field can be seen from (Esh to be 
expressible as 



Boo — {B) 

9 



(52) 



@gi 



ives 



,,,_{B)^ {B){w) (^.)^ , {w)U 
Since, by (p3) the ratio U/N will be given by the formula 



N 



*^o C/.. 



Stt 



N 



(54) 



in which all dependence on Boo and thus also on (B) has 
canceled out, it can immediately be seen that the con- 
ventional definition (which is the same as the definition 
adopted in Ref. pi) 



H = At: 



d{£) 
d{B) 



(55) 



for the effective magnetic field strength H will simply 
give 



H = Bo 



(56) 



i.e. H is directly identifiable with the London field. The 
corresponding magnetic polarization A^, as defined in the 
usual way by 



{B) = H + AttM 
will be expressible as 

(w) _ vNh 



M 



Airq 2q 



(57) 



(58) 



VI. EXTERNAL SOLUTION 

In the previous sections we have been able to establish 
very useful properties concerning the energy density per 
unit length of a vortex, without needing the specification 
of a internal structure. This section and the next one will 
consider this question and show in particular how the un- 
specified parameter of the previous section, the axis value 



of the relative magnetic field, or equivalently the effective 
radius R, can be explicitly computed depending on the 
modelization of the internal structure. We shall focus 
here on the solution outside the core, which is always the 
same up to a normalization constant, that can be deter- 
mined only with the knowledge of the internal structure. 
This will be the purpose of the next section. 

In the region outside the core, i.e. in the range r > ^, 
where the the number density Us is uniform with value 

ris = Uoa , (59) 

the equation obtained from (^IJ) and ( |l4[) by eliminating 
B will have the form 






-^-ft-">'^«' 



where A is a fixed lengthscale given by 



A^ 



m 



A-Kq^Uo 



(60) 



(61) 



and which is called the London penetration length. 

The equation (|6y) is of the well known Bessel-type, 
whose most general asymptotically bounded solution is 
expressible in the form (see Ref. Q) 

V = CKi{x}, (62) 

where the independent variable x is defined by 

x=^-, (63) 

C is a normalization constant and Ki is a modified Bessel 
function. 

It follows immediately from the flux quantization con- 
dition (|l^) that the magnetic potential deviation defined 
by (H) will be given by 



A= CK,{x} + ^— 

q Zirr 



(64) 



where $ is the magnetic fiux integral given by (Eq). 

Using the fact that Ki is related to the Bessel function 
^oby 

Ki = -K , ^0 - -K[ - x'^Ki , (65) 

where a prime stands for differentiating with respect to 
the argument a;, it is straightforward to obtain the corre- 
sponding solution of (113) for the magnetic field deviation, 
which will be expressible in the simple form 



B 



qX 



CKo{x} . 



(66) 



The external configuration for the magnetic vortex has 
thus far been determined up to the normalization con- 
stant C. It will be seen in the next section how, on the 
basis of a suitable ansatz for the radial dependence of 
Ug, the solution inside the core can be used to fix this 
constant C, and thus to determine completely the con- 
figuration of the vortex. 



VII. INTERNAL SOLUTION 

Instead of directly specifying the way in which the 
number density Hs varies from zero on the axis (r — 0) 
to its external value rioo at the core radius (where r = ^, 
i.e. where x — x = ^/A), it is more convenient to work 
with an ansatz based on an explicit prescription for the 
current magnitude j, which will have a qualitatively sim- 
ilar behavior in the core, ranging from zero on the axis to 
a value j at the core radius, that according to (12) and 
( |62| ) will be given by 



j = qUaoCKi 



(67) 



using the obvious abbreviation iiTi = Ki{x}. The current 
in the core, i.e. where x < x, will therefore be expressible 
in the form 



J = crj, 



(68) 



where cr is a dimensionless function of x that is required 
to vanish, a ~ 0, for a; = and to increase to unity, ct = 1 
where x = x {cr plays here the role of the structure func- 
tion mentioned in the introduction). For any suitably 
prescribed function a with these properties, there will be 
corresponding functions, x ^-nd d say, that are defined 
by the requirement that they too should vanish on the 
axis, i.e., x = ^ = Ofora; = and by the requirement 
that they should be obtained in the region < a; < a; as 
solutions of the differential equations 



a; -— = 2ct , 
ax 



X — ; — — ixx ■ 
ax 



(69) 



In terms of such a set of functions, the relevant solution 
of equation ( [ll| ) will evidently be given by 



B = Bo~27rCjX, 



(70) 



and the corresponding solution of (10) will be given by 



A=^XBox-UejC- 



(71) 



The requirement that the magnetic field should be con- 
tinuous (so that the current density j remains finite) en- 
tails that the internal solution ( 7C ) should match the cor- 
responding external solution (36) where a; = a:, so we 
obtain a boundary condition of the form 



Bo = ^C{Ko + ^xK,) , 



(74) 



while the required normalization constant C is finally 
obtained in the form 



C = 



Nh 



XmlC 

in terms of a dimensionless quantity that is given by 

.3 / 4 _ 



(75) 



^=(l + li2^. 



1^2r> 



)xKi + ^x^Ko = 



X 



eKi + -K2 

X 



(76) 



in which the only dependence on the internal structure is 
that embodied in the dimensionless number e which can 
be seen to be given in terms of the boundary values % 
and C, of the functions x ^^'^ C by the simple formula 



2X-C. 



(77) 



This quantity e can be interpreted as the value at the 
core boundary a; = a; of a function e of a; given by 



e = 2^a;-C, 



(78) 



in terms of which the solution for the potential difference 
A will be given by 



mC , ~ x^ ~ , 



(79) 



The corresponding expression for the magnetic field ex- 
cess will have the form 



B^—{lU + -K,{x-x)) 



(80) 



Using the solution ([75|) for C, the central value needed 
for the energy formula (p9) can be seen to be obtainable 
as 



Bo^JS^(2 •*^i(8-C^^) 



qX^x^ 



4/C 



(81) 



VIII. EXPLICIT ENERGY CONTRIBUTIONS 



qX 



(72) 



while the corresponding continuity requirement for the 
potential gives a second boundary condition of the form 



(J) 2rr7 

T^i qi 



(73) 



This pair of boundary equations can be solved to give the 
central magnetic field difference in the form 



Using the axis-field energy formula (^9|) together with 
the solution (|8l|) for the axis-field Sg, we can imme- 
diately obtain the total "dynamical" energy f/dyn = 

C^mag + C/kin = C^mag + t/kin (which is the total energy 
less the condensation energy): 



where 



U, 



dyn 



f/o 



2^0 + xxKi 
2/C 



(82) 



u„ . (m 



V 47rA / 



(83) 



Using the external solution for B and V, one can obtain 
the external magnetic and kinetic energy contributions, 
defined in ( p7|) and (|3l]), in the form 



U, 



cxt 
mag 



~2^ 



{kl 



Kl) 



and 



rrcxt 



niC \ X'^ 



2q 



-K^ki- {kl - k[ 

X 



(84) 



(85) 



where the core structure dependence is contained exclu- 
sively in the constant C. This constant C, using the 
solution (^, can be related to the (core structure inde- 
pendent) constant f/p defined just above by the simple 
relation 



mC 
~2i 



El 

^2' 



(86) 



where, obviously, all the dependence on the core structure 
is contained in the dimensionless term IC. 

As a consequence, the total external energy contribu- 
tion can be rewritten in the simple form 



U" 



xKqK, 
/C2 



Uo, 



(87) 



which can be decomposed in two similar expressions for 
the magnetic and kinetic contributions. The sum of the 
magnetic and kinetic internal contributions can then of 
course be obtained by using the relation [/'"* — [/"^y" — 

r rcxt 



IX. CONCLUSIONS 



Boo has been established here as a precise mathematical 
relation in the framework only of a particularly simple 
model. The problem of generalization to more sophis- 
ticated models, allowing for compressibility, relativistic 
effects and other relevant complications, remains to be 
dealt with in future work. 
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Let us summarize the results of the present work. We 
have first shown that the contributions linearly depen- 
dent on r^oo in the magnetic and kinetic energies cancel 
each other. As a consequence, we find that the effective 
magnetic field strength H is simply the London field. R is 
to be observed that the extra energy density contribution 
arising from the second term in ( pq ) wou ld g ive an extra 
contribution of the form Boo{w) /Anq in (p3|). By includ- 
ing this extra term - overlooking the fact that, according 
to (p3|), it will be canceled by the second term in the ki- 
netic contribution (pQ), which was not taken into account 
- the analysis by MendellEI provided the erroneous con- 
clusion that there would be no polarization, or in other 
words that H should be identified not with the London 
field but simply with the mean induction, meaning the 
replacement of ( pq ) by iJ = (B). 

The identification (|5l) of H, as given by the conven- 
tional definition (pS), with the asymptotic London field 



